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This work examines the critical anisotropy required for the local stability of the coUinear ground 
states of a geometrically-frustrated triangular-lattice ant iferromagnet (TLA). Using a Holstein- 
Primakoff expansion, we calculate the spin- wave frequencies for the 1, 2, 3, 4, and 8-sublattice 
(SL) ground states of a TLA with up to third neighbor interactions. Local stability requires that all 
spin-wave frequencies are real and positive. The 2, 4, and 8-SL phases break up into several regions 
where the critical anisotropy is a different function of the exchange parameters. We find that the 
critical anisotropy is a continuous function everywhere except across the 2-SL/3-SL and 3-SL/4-SL 
phase boundaries, where the 3-SL phase has the higher critical anisotropy. 



Introduction. Geometrically-frustrated systems ex- 
hibit many novel characteristics including non-colhnear 
ground states and multiferroic properties^. One of the 
best realizations of a geometrically-frustrated triangular- 
lattice antiferromagnet (TLA) is CuFe02 , which contains 
stacked hexagonal planes of spin-5/2 Fe'^+ ions. Accom- 
panied by a phase transition from a coUinear 4-sublattice 
(SL) ground state to a non-coUincar phase^i^i^i^, CuFe02 
exhibits multiferroic properties above a critical magnetic 
field or above a critical concentration of non-magnetic 
Al'^+ impurities, which substitute for the Fe'^^ ions^ii. 
Inelastic neutron-scattering cxpcriments^'^i'i^ on CuFe02 
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FIG. 1: (Color online) The 1, 2, 3, 4, and 8-SL phases for 
the ground states of the geometrically-frustrated TLA. The 
solid black lines denote the magnetic unit cell of each phase. 
Up and down spins are designated by red and blue circles, 
respectively. 



have reported a spin-wave (SW) gap of about 0.9 mcV, 
which decreases with Al doping and may vanisbii upon 
the appearance of multiferroic behavior. Similar behav- 
ior is produced in a model TLA as the anisotropy is 
reducedi^ and spin fluctuations about the 4-SL coUinear 
phase become stronger. In this paper, we evaluate the 
critical anisotropies required for the local stability of the 
coUinear magnetic phases in a model TLA with up to 
third nearest neighbors. As shown elsewheroH, the wave- 
vector of the dominant SW instabilities of a coUinear 
phase coincide with the dominant wave- vector of the non- 
colhnear phase that appears with decreasing anisotropy. 
Therefore, an analysis of the critical anisotropies and 
wave-vectors of a frustrated TLA can provide useful in- 
formation about the non-coUinear phases that appear at 
small anisotropy. 

The coUinear ground states of a TLA with strong 
anisotropy were first obtained by Takagi and MakataM, 
who examined an Ising model with interactions up to 
third nearest neighbors. The ground-state phase diagram 
consists of the five phases sketched in Fig. [2 where the 
energies of these five states are given in Table H] Using a 
Holstein-Primakoff (HP) expansion, we have calculated 
the SW frequencies and critical anisotropies for each of 
these phases. 

The Hamiltonian for a TLA is given by 
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where Si is the local moment on site i, Jij is the in- 
teraction between sites i and j, and D is the single- ion 
anisotropy. Employing a HP transformation, the spin 
operators are given by Siz = S — a\ai, Si+ = V^Sai, 
and Si- = V^Sa]. Expanded about the classical limit 
in powers of l/VS, the Hamiltonian can be written as 
H = E + Hi + H2 + ■ ■ ■ ■ The first-order term Hi van- 
ishes when the spin configuration minimizes the energy 
E. The second-order term H2 provides the dynamics of 
non-interacting SWs. Higher-order terms i?n>2 reflect 
the interactions between SWs. They are unimportant 
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at low temperature and for large 1/S. Similar to Tak- 
agi and Makata, we consider nearest neighbor Ji, next- 
nearest neighbor J2, and next-next-nearest-neighbor J3 
exchange interactions, as sketched in Fig. [TJ 

To determine the SW frequencies Wk, we 
solve the equation-of-motion for the vectors 

r (1) (l)t (2) (2)t 1 1 • 1 1 -t- 

Vk = ["^k 1 ^k J , Cj^ , ...J, which may be writ- 
ten in terms of the 2N x 2N matrix M(k) as 
idv-\^/dt = — [gg. Vk] = A/_(k)vk, where N is the 
number of spin sites in the unit cell. The SW frequencies 
are then determined from the condition Pet [A/ (k) — tJw/] 
= 0. 

Two conditions are required for the local stability of 
any magnetic phase: all SW frequencies must be real and 
positive and all SW weights must be positive. The SW 
weights are coefficients of the spin-spin correlation 
function: 



+ (S-S+(t))} = ^M/(^)<S(c.-4^)), 



(2) 



where s denotes a branch of the SW spectrum and is 
defined as the vector pointing from site i to site j. The 

is) 

weig hts W^°' were evaluated within the HP formalism by 
solving the equations-of-motion for coupled spin Green's 
functionsi^ii^. In zero field, the condition that the SW 
weights are positive for all k is equivalent to the condition 
that all SW frequencies are positive. 

We obtained analytic expressions for the SW frequen- 
cies for all phases shown in Fig. [l]with the exception of 
the 8-SL phase, which was solved numerically. Analysis 
of the SW frequencies yields the critical anisotropy Dc 
and the critical wave-vectors k where the SW frequen- 
cies vanish. To simplify the following discussion, the SW 
and anisotropy coefficients are provided in the appendix. 

1 - Sublattice. The 1-SL phase (Fig. [Ija)) is a fcrro- 
magnet with SW frequencies 



2SiD + Ai^ 



(3) 



Since the 1-SL phase is locally stable for any positive 
value of the anisotropy, Dc — 0- The SW intensity W^.^^ 
is constant throughout k for all interactions. 

2 - Sublattice. For the 2-SL phase (shown in Fig. [^b)), 
the SW frequencies are given by 



UJ, 



(2) 
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The SW weights for the 2-SL phase are 
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From Eq. the condition for the local stability of 

a 2-SL phase is - Aj^ > 0. At Dc, Al^ 



^3k- 



TABLE I: Classical Energies and Critical Anisotropies 
for TLA Sublattices 
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This condition is satisfied when Dc — in most of the 
2-SL phase. But approaching the 3, 4, and 8-SL phase 
boundaries, nonzero anisotropy is required for local sta- 
bility. As shown in Fig. ^a), the critical anisotropy is 
continuous across the 4-SL and 8-SL boundaries, but is 
discontinuous across the 3-SL boundary. 

Upon closer examination (Fig. EJc)), we find that Dc 
depends differently on the exchange parameters in the 
three regions designated by Roman numerals. In region 
21 (bounded by J3 = J2/2, J3 = (9J2 - Ji)/12, and 
J3^Ji/{Ji-'2J2)), 



D 
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64J|J2 +48J|Ji 



272J| 

+36J3JIJ1 - 27J| - (2J3 - J2)C3}, 



(6) 



where 



C = v/(2 J3 + 3J2)2 - 8J3 Ji 



(7) 



In region 211 (bounded by J3 = J2/2, J3 
(8J2 - Ji)/9, and J3 = J|/(Ji - 2J2)), 



J2, J3 



(8) 



Finally, in region 2III (bounded by J3 — J2/2, J3 — Ji/4, 
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FIG. 2: (Color online) SW frequencies at the critical 
anisotropy for 2 (a), 3 (b), 4 (c), and 8-SL (d) phases (in- 
teraction parameters given in the text). All SW instabilities 
occur for kya = except in regions 2III and 81, where they 
occur for kya = O.ISGtt and 0.3827r, respectively. 



and Js = (Ji - J2)/4), 



^(2111) ^ (4/3 + J2- Ji)^ 

2( J2 + 2,h) 



(9) 



The region with no critical anisotropy D^c^^^ is bounded 
by J2 = 0, J3 = (8J2 - Ji)/9, and J3 = (Ji - J2)/4 as 
shown in Figs. [3]Jb) and (c). 

The critical wave- vectors, k for the SW instabilities in 
region 21 are: 



(21,1 



2 arccos 



3 J2 - 2 J3 - C 
8J3 



(10) 



Two other instabihties k^^^^*^) and k*^^''^) are related to 
]^(2i,a) ±7r/3 rotations and can be considered "twins" 

of the k^^^'*^-* instabilities. All three instabilities occur at 

(21) 

the same critical anisotropy Dc ■ For regions 211 and 
III, the SW instabilities occur at 

TT ± 7r/3, 



(11) 



0, 



and 



2 r J2 + Ji 

arccos < — 

V3 I 2( J2 + 2 Ja) 



(12) 
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FIG. 3: (Color online) Critical anisotropics for the TLA 
ground states. The spacing between contours is 0.2 (a), 0.1 
(b), and 0.05 (c). Numbers designate the stable phase and 
Roman numerals designate regions where the behavior of the 
critical anisotropy is distinct. Solid lines denote boundaries 
between phases and dashed lines denote boundaries between 
regions. is continuous across each phase boundary except 
for the 2-SL/3-SL and the 3-SL/4-SL phase boundaries; in 
both cases, Dc is higher for the 3-SL phase. 



along the and ky axis, respectively. 
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Notice find that oi^^ = along the 3-SL/l-SL bound- 
ary. Again, Dc is discontinuous along the 2-SL/3-SL and 
3-SL/4-SL boundaries: the anisotropy required for the 
local stability of the 3-SL phase is three times the crit- 
ical anisotropy of the 2 or 4-SL phases. As discussed 
further below, the discontinuities at the 2-SL/3-SL and 
3-SL/4-SL phase boundaries are related to the distinction 
between the conditions for global and local stability. 

In Fig. El^b), we plot a SW dispersion in the 3-SL phase 
with interaction parameters J2/IJ1I = 0.5, J3/IJ1I = 
—0.5, and Dc/\Ji\ = 2.25. Since the 3-SL phase has a 
net moment, the SW frequencies are quadratic functions 
of k near the instability wave- vectors. 

4 - Sublattice. The SW frequencies for the 4-SL phase 
(shown in Fig. [ijd)) were evaluated in Ref. 
given by 



12l | and are 



FIG. 4; (Color online) SW frequencies and SF intensities 
for the 2, 4, and 8-SL phases at the critical anisotropy with 
■h/\Ji\ = -0.5, J3/IJ1I = -0.25, and D/\Ji\ = 0.125. 



^7k 



+4|A6k^2k - A7ki^2kl 



^2k ^2k 



l/2\ 1/2 



(16) 



IV were evaluated with kyU 
evaluated at kyU 



Figure [^l^a) shows three representative SWs for all 2- 
SL regions. The interaction parameters for region 21 are 
J2/IJ1I = -0.25, J3/IJ1I = -0.12, and DJ\Ji\ = 0.04. 
For region 211, J2/IJ1I = -0.10, J3/IJ1I = -0.05, and 
DJ\Ji\ = 0.325. For region 2III, J2/IJ1I = -0.75, 
J3/IJ1I = -0.125, and Dc/\Ji\ = 0.031. Finally, for re- 
gion IV, the interaction parameters are J2/IJ1I = —1.0, 
J3/IJ1I = -0.125, and Dc/\Ji\ = 0.0. Regions I, II and 

0, while region III was 
0.1867r as explained above. 
In Figs. ISI^b) and (c), we examine the critical 
anisotropy of the 2-SL along the J3/IJ1I = axis. The 
critical anisotropy vanishes for —1 < J2/\Ji\ < —1/8 
but is nonzero outside this region. Therefore, non- 
coUincar phases should appear for J2/\Ji\ < —1 and 
^2/1^1! > —1/8 when D < Dc- This agrees with Joli- 
coeur et ali^, who studied a TLA with nearest and next- 
nearest neighbor exchange interactions and D = 0. They 
obtain a Neel state up to J2/I = —1/8 and an incom- 
mensurate spiral for J2/\Ji\ < -1- Similar results have 
been obtained on square lattice a^^i^^ . 

3 - Sublattice. For the 3-SL phase (shown in Fig. [Uc)), 
the SW frequencies are 



,(3) 



/^cos(6'/3 + 277i/37r) + i?2k/3, (13) 



where m is an integer (0,1,2) distinguishing the three 
separate SW dispersion relations and 



2i?2\ 



9i?2k^3k — 27i?4k I 



1458i? 



3/2 

Ik 



(14) 



The critical anisotropy of the 3-SL phase is independent 
of J2 and given by 



(15) 



The SW weights of the 4-SL phase are 

Wi"^^ = [i?5k(A7k - ^6k) + (i^2k + i^2*k)(^6k - ^7k)' 



iF2^.~F*^)^F2^. + R 



2k 



A. 



6k 



A 



7k, 



^5kV^6k - ^7k - ^5k 

(17) 

As for the 2-SL phase, the critical anisotropy for the 
4-SL phase depends differently on the interaction param- 
eters in two regions, again denoted by Roman numerals I 
and II (Fig. [DJb)). In region 41 (bounded by J3 = J2/2, 
J2 = Ji/2, and J3 = J|/( Ji - 2J2)), 
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(41) 



(4J3) = 



16J| 



64Jf J2 



48 J| Ji 



48JU2J1 



(18) 



+36J3JIJ1 



27J| 



(2 J3 - J: 



and in region 411 (bounded by J3 

J3 = Jl/{Jl-2J2)), 



J2/2, J2 = 0, and 



i5(4n)^2J2-^J3-^Ji. 



(19) 



The critical wave- vectors for the 4-SL phase are the same 
as those in the respective region of the 2-SL phase, includ- 
ing the multiple instabilities in region 21: k^"'^''*^ = 
k(4i,b) ^ i^(2i,b)^ ^j^j 1^(411) ^ ]^(2ii)^ Figure ll^c) shows 

two representative SWs for regions 41 and 411 with kya = 
0. The interactions parameters for region 41 are J2/ 1 >/i | = 
-0.439, J3/IJ1I = -0.570, and Dj\Ji\ = 0.105. For re- 
gion 411, they are J2/I Ji| = -0.25, J3/I Ji| = -0.5, and 
Dj\Ji\ = 0.25. 
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FIG. 5: (Color online) (a) Location of SW instability A — a\kx — 7r| along kya = in region 41 for fixed values of J2/jJi|. As 
IJsl/IJil increases along J2/IJ1I = —0.5, A asymptotically approaches tt/3. (b) Plot of A in region 41 along the 2-SL/4-SL 
boundary J3/J2 ~ 1/2. The cusp in A occurs at J2/IJ1I = —1/3 where the SW instability occurs at tt. 



Figure |4] shows the SW frequencies and SF intensi- 
ties for the 2, 4, and 8-SL phases at the triple point of 
the phase diagram. The intensities for the 8-SL phase 
were determined numerically. As shown in Fig. [4l each 
phase becomes unstable at D/\Ji\ — 0.125, where the 
SW intensity for each phase peaks at the same wave- 
vector. This wave-vector corresponds to the ordering 
wave- vector of the non-collinear phased that appears at 
small anisotropy. Because the 2, 4, and 8-SL phases all 
have zero net moment, their SW frequencies are linear 
functions of k around the wave-vectors of the instabili- 
ties. 

The 4-SL phase is of particular interest since it is the 
known ground stated of CuFe02. Fits of the experi- 
mental SW frequencies^iii^ of CuFe02 have determined 
the ratios of exchange parameters Jil\J\\ ~ —0.44 and 
>^3/|>^i| ~ —0.57, which lies within region 4L Conse- 
quently, we have studied the SW frequencies of the 4- 
SL phase more closely. Figure shows the behavior 
of fci'*^'''' along various cuts through region 41 of phase 
space. Since the SW frequencies are symmetric about 
the midpoint of the Brillouin zone ott, we consider the 
quantity A = ajfcj^- 7r|. As J3/I Ji| increases in region 41, 
A asymptotically approaches 7r/3, which is the constant 
value of A in region 4IL For small values of J3/I Ji|, the 
wave- vector instabilities approach tt as J2l\J\\ increases, 
equal TT for J2/I Ji| = —1/3, and then move away from tt 
as J2l\J\\ approaches zero; this behavior is shown along 
the 2-SL/4-SL boundary in Fig. [^b). 

8 - Sublattice. For the 8-SL phase (shown in Fig. [He)), 
we have determined SW dispersion relations numerically. 
The critical anisotropy values for this phase are shown in 
Fig. [nja). Notice that Dc has a cusp dividing the phase 
into regions 81 and 811 (Fig. [2Ib)), separated by J3 = 
J2/2. Looking more closely at the numerical results, the 
critical anisotropics in the 8-SL regions arc closely related 



to those of their respective neighbors and are given by 
=i?f")+4J3- Ji, (20) 

i^(8ii) ^^(4I)+2J2- Ji, (21) 

which clearly show that the critical anisotropics are con- 
tinuous across the phase boundaries. In region 811, the 
wave-vector instabilities occur for ky = (as in region 
41); in region 81, the wave- vector instabilities occur for 
non-zero ky (as in region 2III). Figure^d) shows two rep- 
resentative SWs for regions 81 and 8n. The interactions 
parameters for region 81 are J2/IJ1I = —1.5, J3/IJ1I = 
-0.50, and Dc/\J\ = 0.25. For region 811, they are 
J2/IJ1I -0.75, J3/IJ1I = -0.50, and Dj\Ji\ = 0.62. 
Whereas kyU — for region 811, kya = 0.3827r for region 
81 as explained above. 

To better understand the discontinuities along the 2- 
SL/3-SL and 3-SL/4-SL phase boundaries, we consider 
the relationship between local and global stability. Our 
SW calculations only guarantee the local stability of each 
collinear phase. But even when a phase is locally sta- 
ble, it can still be globally unstable to a lower-energy 
spin configuration. Hence, the critical anisotropy Dc 
for global stability must be greater than or equal to the 
critical anisotropy Dc for local stability. Unlike Dc, Dc 
must also be a continuous function of Ji, J2, and J3. 
So when Dc is discontinuous, the phase with the lower 
critical anisotropy cannot be globally stable. Since the 
3-SL has a higher critical anisotropy along the 2 and 4-SL 

boundaries, the 2-SL and 4-SL phases cannot be globally 

(211) (3) 
stable along those boundaries when Dc < D < Dc 

or ZJi^^^' < D < D^\ Therefore, our results for the local 
stability of the collinear phases also has implications for 
the global stability of those phases. 

Conclusion. We have examined the critical anisotropy 
for a geometrically-frustrated TLA. Based on the Takagi- 
Makata phase diagram, wc calculated the SW frequencies 



for all five phases. Imposing the two conditions for lo- 
cal stability, we obtained the critical anisotropies and 
wave-vector instabilities for all phases as functions of 
the exchange interactions. Surprisingly, these results are 
highly dependent on the longer-range exchange interac- 
tions and most phases break into several regions where 
the anisotropy has a distinct dependence on the exchange 
parameters. As discussed for the 2-SL and 4-SL phases, 
the critical anisotropies and wave-vectors for the local 
stability of the coUinear phases provides useful informa- 
tion about the non-coUinear phases that appear at small 
anisotropy. We have also shown that the discontinuity of 
the critical anisotropy at the 2-SL/3-SL and 3-SL/4-SL 
phase boundaries has implications for the global stabil- 
ity of the 2-SL and 4-SL phases with the smaller critical 
anisotropies. 
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^3k = -^1 ( cos (k • d2) + cos (k • da) 
+J2 ( cos (k • + cos (k • dg 



The 3-SL phase coefficients are 

^Ik = ^2k ^ 3i?3k;, 
i?2k = 2A4k -I- Ask, 

-Rsk = ^Ik + 2A4kA5k + l-Fikp, 



(A3) 



(A4) 



(A5) 



(A6) 



-R4k = (Ask - 2A4k)|^^ikP - A^^Ask - - F*^, 



= 2D + 2J2(3 - cos(k • d4) 
— cos(k • d.s) — cos(k • dg)). 



(A7) 
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APPENDIX A: SPIN- WAVE AND ANISOTROPY 
COEFFICIENTS 

This Appendix provides the coefficients that enter the 
SW frequencies and weights for each phase. The coeffi- 
cients for the 1-SL or ferromagnetic phase are 



Aik = 3(Ji + .h + J3) 



- Ji ^ cos (k • di) + cos (k • d2) + cos (k • 



— J2 ( cos (k • d4) + cos (k • ds) + cos (k • dg) 



- J3 cos (2k • di) + cos (2k • d2) + cos (2k • dg) j , 

(Al) 

where di = ax, d2 = 1/2 ax+ V3/2 ay , da = —1/2 ax-|- 
-\/3/2ay, d4 = 3/2 ax + \/3/2 ay , ds = ay , and 
dg = -3/2ax + \^/2 ay. 

The 2-SL phase coefficients are 

A2i.^D + 3 J3 

-Ji('cos(k • di) + 1) - J2(cos(k • ds) -f 1 



- J3 ( cos (2k • di) + cos (2k • d2) cos (2k • dg) j , 

(A2) 



^k = 6 Ji + 6 J3 - Asu, 



i^ik = Ji(e 
+J3{e 



22k-d2 I ^ — 2zk-di i ^ — 2ik-d3 



As in Ref . (fl^ ) . the 4-SL phase coefficients are 

^5k = (-F2k + ^2k ~ 2(^2*^ -f- 2A6kA7k)^fk 



(A9) 



(AlO) 



+4(AL + ^?k)l^2kP - 4AgkA7kF2*2 



Aei, = D- Ji + J2(l - cos(k • ds) 



1/2 



-J3(l + cos(2k-di)), 



(All) 



(A12) 



^rk = — cos(k • di)(ji -f 2J3Cos(A/3k-d5)), (A13) 

^2k = - C0s(k • d5/2) _^ J2g-3ik.di/2^ ^ 

(A14) 
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